Astronomy & Astrophysics manuscript no. astroph © ESO 2008 

October 14, 2008 



Extrasolar planets and brown dwarfs around A-F type stars * 

VI. High precision RV survey of early type dwarfs with HARPS 

A.-M. Lagrange 1 , M. Desort 1 , F. Galland 1 , S. Udry 2 , and M. Mayor 2 

1 Laboratoire d'Astrophysique de l'Observatoire de Grenoble, Universite Joseph Fourier, BP 53, 38041 Grenoble, France 

2 Observatoire de Geneve, 51 Ch. des Maillettes, 1290 Sauverny, Switzerland 

Received date / Accepted date 

ABSTRACT 

Aims. Systematic surveys to search for exoplanets have been mostly dedicated to solar-type stars sofar. We developed in 2004 a 
method to extend such searches to earlier A-F type dwarfs and started spectroscopic surveys to search for planets and quantify the 
detection limit achievable when taking into account the stars properties (Spectral Type, v sin i) and their actual levels of intrinsic 
variations. We give here the first results of our southern survey with HARPS. 

Methods. We observed 185 A-F (6 - V in the range [—0.1; 0.6]) stars with HARPS and analysed them with our dedicated software. 
We use several criteria to probe different origins for the radial-velocity variations - stellar activity (spots, pulsations) or companions: 
bisector shape, radial-velocity variations amplitudes and timescales. 

Results. 1) 64% of the 170 stars with enough data points are found to be variable. 20 are found to be binaries or candidate bina- 
ries (with stars or brown dwarfs). More than 80% or the latest type stars (once binaries are removed) are intrinsically variable at a 
I 2 m s~'precision level. Stars with earlier spectral type (B -V < 0.2) are either variable or associated to levels of uncertainties com- 

parable to the RV rms observed on variable stars of same B - V. 2) We have detected one long-period planetary system (presented 
i | in another paper) around an F6IV-V star. 3) We have quantified the jitter due to stellar activity and we show that taking into account 

this jitter in addition to the stellar parameters (spectral type, v sin f), it is still possible to detect planets with HARPS with periods of 3 
O | days (resp. 10 days and 100 days) on 91 % (resp. 83 %, 61 %) of them. We show that even the earliest spectral type stars are accessible 

to this type of search! , provided they have a low projected rotational velocity and low levels of activity. 4) Taking into account the 
present data, we compute the actually achieved detection limits for 107 targets and discuss the limits as a function of B - V. Given 
the data at hand, our survey is sensitive to short-period (few days) planets and to longer ones (100 days) at a lower extent (latest type 
stars). We derive first constrains on the presence of planets around A-F stars for these ranges of periods. 
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V^jQ 1 . Introduction terns as well as the lack of solid material close to the star could 

£C) reduce the number of planets. Clearly, several parameters proba- 

\0 Slnce the dlscover y of the first exoplanet around a solar-like bly impact the occuren ce and properties of planets around mas- 

^j- star in 1995, more than 250 planets have been found by radial- siye ^ and they haye not t been fully explored 
velocity (RV) surveys (Jean Schneider, http://exoplanet.eu). 

These surveys have generally focused on late-type stars (later The data to test the models are still quite limited as the 

5rg than F8). However, knowing about the presence of planets or largest and earliest, now long-lasting surveys had focused on 

q brown dwarfs (hereafter BDs) around more massive objects is solar type, main-sequence (MS) stars. In recent years, some ef- 

• • mandatory if one wishes to investigate the impact of the mass of forts have been made nevertheless to search for planets around 

the central stars on the planetary formation and evolution pro- stars with various masses: less massive, M-type stars on the one 

cesses, hand, and more massive stars on the other hand. The search 

There are theoretical indications that the mass of the planets for P lanets around M stars seems sofar to confirm the pre- 

c3 incre ases with the mas s of the parent star, at least for low mass viou sly mentionned expectations from the oretical works (see 

stars (llda& Lin 2005b and that the frequency of giant planets e -8-> IBonfils et al. 2005J |Butler et al. 2006). As far as massive 

increases linear ly with the parent star ma ss for stars between stars are concerned, the available observations are still veiy 

0.4 and 3M ( |Kennedy & Kenyon 2008| l, with e.g., 6% fre- limited. Massive MS stars have been removed from early sur- 



quency of giant planets around 1 M and 10 % frequency around ve y s as it was generally thought that their spectra (few lines, 
1.5 M . More numerous and massive planets are consistent with usually broadened by stellar rotation) would not allow planet 
what we could expect from a disk surface density increasing with detection and indeed, the classical RV measurements technics 
stellar mass. On the other hand the shorter lifetimes of the sys- ( ba sed on the cross correlation of the actual spectra with a bi- 
nary spectral mask con-esponding to a star with an appropriate 
Send offprint requests to: A.-M. Lagrange, spectral type and vsin/ = km s" 1 ) fail to measure the RV of 
e-mail: anne-marie.lagrange@obs.ujf-grenoble.fr these stars. This lead some groups to study instead "retired" 
* Based on observations collected at the European Southern early-type, either low-mass (< 1.6 M ) giants, intermediate- 
Observatory, Chile, ESO 075.C-0689, 076.C-0279, 077.C-0295, 078.C- mas s (1.6-2M ) subgiants or clump giants (1.7-3.9 M ) (see 
0209, 080.C-0664, 080.C-0712. e.g. JHatzes 20 05 Niedziels kTet al. 20071 fjohnson et al. 2006t 
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Johnson et al. 20071 |Lovis and Mayor 2007| ISato et al. 2 008 ). 
These stars indeed have cooled out and also rotate more slowly 
due to coupling of stellar winds and magnetic fields; they there- 
fore exhibit more numerous, narrower lines, which is adequate 
for classical RV measurements technics, and their level of ac- 
tivity (jitter) is relatively low (10-20 ms _I ; Hekker et al. 2006) 
for giants and 10 ms" 1 for subgiants ( Joh nson et al. 20071 and 
ref. therein; ISato et aL~2 008 ). The data available today are still 
limited compared to those available on solar-type stars and less 
than 20 planets have been found sofar in total around these 
evolved stars. So far, the planets found around K subgiants stars 
with M > 1.5Mq are located at distances larger than 0.8 AU 
dJohnson et al. 200Tb : this led these authors to conclude that 
close-in planets are rare, in agreement with some theoretical pre- 
dictions on disks depletion timescales (Burkert and Ida 2 0071 1. 
However, the impact of the post MS evolution of the stars 
on closer-in planets has not been explored yet for these stars. 
Concerning giant stars, all planets found so far have relatively 
long periods, the closest ones beeing reported at less than 0.7 AU 
from 2 giants, in addition to a previously reported planet at 
0.7 AU by Sato et a l. 20031 Numerical simulations by the same 
authors suggest that planets with orbits inside 0.5-1 AU around 
2-3 Mq stars could be engulfed by the central stars at the tip 
of RGB due to tidal torque from the central stars. According to 
them, if one assumes then that most of the clump giants are post 
RGB stars, there is then a risk that closer planets, if present be- 
fore had disappeared since the star evolved. In summary, even- 
though there are some hints that hot Jupiters are not present 
around retired stars, it is recognized that data are still needed 
to definitely confirm this point. Also, as stellar evolutionary pro- 
cesses may have impacted the presence of planets close to the 
stars, it is acknowledged that data are needed on A-F MS stars 
(see eg IBurkert and Ida 20071 see also |Li, Lin and Liu 2 008 ). 
We note finally that short-period planets have indeed been 
found around F5-F6 MS stars through transit observations (see 
http://exoplanet.eu). 

We developed a few years ago a dedicated software to ex- 
tract the RV data around early type MS stars. The method con- 
sists in correlating, in the Fourier space, each spectrum and a 
reference spectrum built by summing-up all the available spec- 
tra for this star. We had shown earlier that with this approach, 
and taking into account the stars B - V and projected rotational 
velocities, it is possible to find planets around A-F type stars 
dGalland et al. 2005al ). However, the price to pay is that more 
measurements are needed to find planets around A and early 
F dwarfs than for late F and G-K dwarfs, because of the rela- 
tively higher uncertainties in the RV measurements due to higher 
v sin i and higher effective temperature, and also because of the 
possible presence of pulsations or spots in the case of late F 
stars, the impact of which has not been quantified so far. (Note 
that spots or pulsations become also a limiting factor in the 
case of later-type stars, as well as pulsations, if one looks for 
low mass planets.) We started then systematic searches for low- 
mass companions to A-F type stars, with HARPS in the south- 
ern hemisphere and with ELODIE and then SOPHIE at OHP in 
the northern hemisphere. We found so far a 9.1 M Jllp (minimum 
mass) planet orbiting (a — 1.1 AU) an F6 V star, with v sin i = 
12kms _1 (Gal land et al. 2005bl . Very promisingly, we also de- 
tected a 21 M Jup brown dwarf orbiting (a = 0.2 AU) a pulsating 
A9V star with vsin/ = 50kms _1 dGalland et al. 20061 ): notice- 
ably, in that case, we could disentangle stellar pulsations from 
the presence of a low-mass companion. 

In parallel, we developed detailed simulations of stellar ac- 
tivity (spots) in order to estimate more quantitatively than what 



was available so far dSaar & Donahue 19971 [Hatzes 2002} the 
impact of such stellar activity on RV data and other observ- 
ables (bisectors, bisectors velocity-span, photometry). We have 
showed that if the star vsin / is smaller than the spectrograph 
spectral resolution, depending on their location with respect to 
line of sight, depending also on their size, spots with realis- 
tic sizes can produce RV variations and bisector velocity-span 
variations quite similar to those of low mass planets. Hence, 
low amplitude (level of typically 20ms~' or less) planetary- 
like RV and bisector velocity span variations cannot alone defi- 
nitely prove the presence of planets around low v sin i G-K stars 
dDesort et al. 20 07 ) and additional criteria are mandatory to rule 
out spots: photometry, activity evaluation down to levels relevent 
to explain the amplitude of RV variations, precise knowledge of 
the star rotaional period, etc. The situation in that respect is much 
more favorable in the case of earlier-type stars as they rotate sta- 
tistically faster and hence the bisector criteria can apply. 

The present paper is devoted to our southern hemisphere sur- 
vey. The sample, observations, measurements and diagnostics 
are provided in Sect. [2] The results concerning the stellar vari- 
ability as well as the quantitative impact on planet detectability 
around the early-type stars are presented in Sect. [3] Finally we 
give and discuss in Sect. |4] the detection limits obtained in the 
present survey. 



2. Sample, observations and measurements 

2.1. Sample 

Our HARPS sample is limited to B8 to F7 dwarfs. The limit 
in spectral type (ST) at F7 is set because the surveys us- 
ing masking technique generally start with stars with ST later 
than F8. The limit at B8 is set by the precision that can be 
obtained with our method on stars given their ST and their 
v sin i ( Gallan d" et al. 2005al ): the detection limit of stars with ST 
earlier than B8 does not fall into the planet domain. Our survey is 
also volume-limited with an upper limit at 67 pc for the B8-A9 
dwarfs and at 33 pc for the F0-F7 dwarfs. The distance was taken 
from Hipparcos catalogue, and stars with distance uncertainties 
larger than 20 % were removed. The difference in distance for 
both spectral types comes from the fact that we wanted to have 
roughly the same number of A and F stars. The dwarf nature was 
selected by selecting stars with absolute magnitudes within 2.5 
magnitude from the Main Sequence. 

Spectroscopic binaries as well as close visual bina- 
ries with separations smaller than 5" known at the be- 
gining of the survey from Coravel or Hipparcos data 
were removed. Confirmed 5 Scuti (from Rodriguez et al. 2000 ) 
and y Doradus type stars (from Mathias et al. 2004 and 
http://astro.univie.ac.at/dsn/gerald/gdorlist.html) were also re- 
moved as they are known to produce RV variations over hours 
to a few days periods due to pulsations. Finally we removed as 
well Ap and Am stars, which present spectral anomalies and are 
often associated to binary systems. This removes a number of 
late A - early F type stars, crossing the 5 Scuti and y Doradus in- 
stability strip. We ended up with 207 stars with ST between B8 
and F7, and B — V ranging respectively between -0.1 and 0.58, 
corresponding to mass ranging between 1.3 and 3.5 Mq. Note 
that we have a relatively smaller number of stars in the [0.2; 0.4] 
B - V range (i.e., roughly, between 1.8 and 1.4M©) as we re- 
moved the known 8 Scuti and y Doradus stars). 
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2.2. Observations 

185 stars have been observed between August 2005 and January 
2008. Figure[T|shows their position in the HR diagram. It can be 
seen that our survey fills a domain of the HR diagram that was 
not covered yet. 

We usually recorded 2 consecutive high resolution (R 
1 15000) spectra each time we pointed to the star (each point- 
ing is hereafter refered to as one epoch). The spectra cover a 
wavelength range between 3800 and 6900 A. As far as possible 
we moreover tried for a given object to record data at two or 
three different times during one night, in order to identify possi- 
ble high-frequency RV variations. We also tried whenever pos- 
sible to record data on two or three consecutive nights. The time 
baseline for a given star varies between 5 days and more than 
800 days. Only a few (11) stars have been observed during one 
night only, but for 15 stars we got only 4 good quality (i.e., with 
an absorptiorQsmaller than 2) spectra or less. We ended up with 
170 dwarfs for which we recorded 6 or more good quality spec- 
tra. 45 have B - V < 0. 1 ; 72 have B-V between 0. 1 and 0.4 and 
53 have B - V > 0.4. We will restrain hereafter our study to those 
stars. 

Typical exposure times ranged between 30 s and 15min de- 
pending on the star magnitude and on the atmospheric condi- 
tions. 

Table[T| provides the 170 targets observed, together with a 
number of relevent information on the stars (ST, v sin i, B - V) 
and on the data obtained as well as various measurements (see 
below). 
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Fig. 1. Observed stars in a HR diagram. We also plotted the 
dwarfs and/or (sub-) giants surveyed either with the Coralie 
spectrograph or by Johnson et al, 2006. Our targets cover a do- 
main that was not surveyed yet. Note the relative lack of objects 
in the [0.2; 0.4] B-V region, due to selection effects (see text) 



1 Magnitude difference between the observation and the 
one that gives the best signal-to-noise ratio (SN): abs(i) = 
61og 10 [6 rv (i)/ min(e lv )], where e„(l) stands for the uncertainty as- 
sociated to the measurement of the observation (i) for the considered 
object, and min{e„) is the smallest value of uncertainty obtained for 
this object. The measured uncertainties (cf Galland et al. 2005aj take 
into account the photon noise + instrumental uncertainties. 



2.3. Measurements 

2.3.1. Radial velocities 

The extraction of the radial velocities is fully described in 
IGalland et al. 2005al Briefly, for each star, we build a first esti- 
mate of the reference spectrum which is the average of the spec- 
tra recorded and reduced via the STS HARPS pipeline. We then 
compute a first estimate of the RV for each spectrum, by corre- 
lating in the Fourier space each spectrum and this first estimate 
of the reference spectrum. We then build a final reference spec- 
trum by averaging the spectra once shifted from their measured 
RV. For each spectrum we finally measure the RV velocity with 
respect to this reference spectrum. We also measure the uncer- 
tainties associated to each RV measurement. 

Note that to build up the reference we compute the x 2 of 
each spectrum compared to the first estimate of the reference 
spectrum. Most of the time, the x 2 found is much less than 10. 
Whenever a higher x 2 was found, we checked the spectra. In 
such cases, either they were due to bad observing conditions or 
technical problems and were not kept to build the reference spec- 
trum (this actually happened quite rarely as we already selected 
spectra with acceptable absorptions) or they were associated to 
lines deformations indicative of a type-2 binary. 

2.3.2. CCFs, bisectors and bisectors velocity-span 

Whenever possible (see below) we computed for each target 
the resulting cross-correlation functions (CCFs) and the bisec- 
tors velocity-span (see for their definition Galla nd et al. 20 05a). 
Indeed, the bisector and bisectors velocity-span are very good di- 
agnostics of stellar activity (spots, pulsations) provided 1 ) they 
can be measured (see below), and 2) the star projected rota- 
tional velocity is larger than the instrumental resolution (see 
IDesort et al. 20071 

The uncertainty associated to the bisectors velocity-span de- 
pends directly on the projected rotational-velocity and/or their 
spectra type. Indeed, the number of lines used to compute the 
CCF depends on these two parameters (much more than on 
the signal-to-noise ratio). For stars with high vsku (typ. > 
150 km s~') and/or B - V < 0.1, the number of lines may be 
quite low (30-50) whereas for late-type stars with moderate 
v sin i (10-20 km s~'), the number of lines used is a few hundreds 
(up to about 1 000). When the bisectors were computed, we then 
attributed quality flags to the bisectors velocity-span measure- 
ments, respectively: Good, Acceptable, Bad, corresponding to 
numbers of lines respectively > 100, 40-100, and < 40. 

2.4. Diagnostics for the classification of variable stars 

Variable stars are defined as having a RV standard deviation 
(rms) larger than twice the RV uncertainties and a total RV am- 
plitude larger than 6 times the RV uncertainties. RV variations 
can a priori be due to the presence of a companion (star, brown 
dwarf, planet) or to intrinsic variations of the star (spots, pulsa- 
tions). It can also be a combination of those different origins. 

2.4.1. Binarity 

We first checked those stars with high^- 2 (> 10) and looked for 
line deformations indicative of spectroscopic binaries. Figure [2] 
provides an example of a binary SB2 identified on the basis of 
the^ 2 , HD2885 (A2V; vsim = 40kms^). It has to be noted 
that in such cases, the RV values measured are not any longer 
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valid, as our RV extraction method assumes that all lines in a 
given spectrum originate from the same object. 

For the rest of the variable stars, we tried to identify binary 
stars among the stars for which the RV amplitude can be ex- 
plained by the presence of a stellar or BD companion. To do 
so, using rough estimations of the star masses via their B — V, 
we computed the RV amplitude 2x/f2d expected from the pres- 
ence of a 13 Mj up body orbiting with a period of 2 days and the 
RV amplitude 2x A^ood expected in the case of a 200-days pe- 
riod and compared these quantities to the observed RV ampli- 
tudes, once corrected from the RV variations observed within a 
night (in practice, over a few hours), as the variations occuring 
within a few hours are assumed to be due to stellar origin, see 
below. Quantitatively we define as "in-night" RV amplitude for 
a given object, the amplitude of the nightly RV variations. We 
computed then the following quantities: 7?2d = (observed RV am- 
plitude - "in-night" RV amplitude)/ 2 x^a and /?200d = (ob- 
served RV amplitude - "in-night" RV amplitude)/ 2x/f2ood, *° 
be used as thresholds to identify the binaries. Note that we chose 
2 and 200-day periods as they are quite relevent given our tem- 
poral sampling and our average time baseline. 

For those variable stars for which we could compute a CCF 
and test the relation between bisectors velocity-span and RV 
variations, we selected those that show a simple, flat bisectors 
velocity-span, i.e., values of bisectors velocity-span arranged 
horizontally in a (RV; bisectors velocity-span) diagram (this cor- 
responds to stars for which the ratio of the amplitude of the bi- 
sector velocity span to the RV amplitude is smaller than 0.2) 
and for which Rjd > 2 or 7?2ood - 2; we regard them as unam- 
biguous binaries. Figure [3] provides an example of a variable 
star (HD 68456; F5V; vsin i = 12kms _1 ) for which the bisectors 
velocity-span clearly indicates the presence of a companion and 
the observed RV amplitude once corrected from in-night varia- 
tions can be due to a ^ 0.1 M0 stellar companion (see below). 
Note that this star was also recently classified as a binary on the 
basis of astrometric data (see below; Goldin & Maka rov 20071) . 

Some variable stars show a bisectors velocity-span which is 
either partly flat and partly vertical or partly flat and partly in- 
clined, indicating that most probably they are binaries and at the 
same time pulsating or active (see below). We classify those ob- 
jects with 7?2d or ^200d > 2 as strong binaries candidates. An 
example, HD 19545 (A3V; vsin; = 80km s _1 ), is provided in 
Fig. [4] Figure|5]shows for comparison the case of a pulsating star 
for which we artificially simulated an additional companion star. 
The generated RV and span curves of the pulsating star and pul- 
sating star plus stellar companion are comparable to those found 
in the case of HD 19545. Note that we do not have quantitative 
criteria to identify those "composite bisectors velocity spans"; 
this is why we classify the candidates as strong candidates rather 
than unambiguous binaries. 

Finally, for the rest of the stars, we flagged as binary can- 
didates those stars with 7?2d or ^2ood ^ 4. Note that for these 
stars we conservatively adopted a more stringent threshold for 
/?2d or 7?2ood as we lack of additional indication of the presence 
of companions, and we know that these stars may be intrinsi- 
cally variable; we thus took into account the fact that the actual 
amplitude RV variations due to the pulsations may be larger than 
the one measured on our set of data. This ensures that most of 
the observed RV amplitudes is due to a perturbation by a BD or 
a star. Figure|6]gives an example of such a star. 




Radial Velocity [km s '] [km s ] ] 



Fig. 2. Example of an SB2 binary, HD2885 (A2V; vsin/ = 
40kms _1 ). RV curve (upper right), CCFs (lower right), bisec- 
tors (upper left), and bisectors velocity-span (lower left). The 
CCF is clearly variable and indicative of an SB2 binary. 




Fig. 3. Example of a binary identified by a flat bisectors velocity- 
span diagram and high amplitude RV variations , HD 68456 
(F5V; vsin / = 12kms _I ). RV curve (upper right), CCFs (lower 
right), bisectors (upper left), and bisectors velocity-span (lower 
left). The mass of the companion falls in the stellar domain (see 
text). 



2.4.2. Planets 

Those stars that shows at the same time signs of RV variability 
with low amplitudes and a flat bisectors velocity-span diagram 
are very good candidates for hosting planets. Note that in some 
cases, stars showing composite bisectors velocity-spans diagram 
with /?2d or /?200d larger than 2 can still be "good" candidates 
for hosting planets. In such cases, the total RV amplitude is not 
dominated by the planetary signatures but but stellar variability 
{e.g., spots). 

2.4.3. Intrinsically variable stars 

In the case of spots, and provided the star v sin i is higher than 
the instrumental resolution, the bisector shape is very peculiar 
and the bisectors velocity-span variations are correlated to the 
RV ones (see Des ort et al. 200 7). In a (RV; bisectors velocity- 
span) diagram, the bisectors velocity-span values are arranged 
either as an inclined "8" shape, or along an inclined line (so 
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Fig. 4. Example of a star whose RV variations are most probably 
due to both pulsations and binarity, HD 19545 (A3V; vsin/ = 
80km s~'). RV curve (upper right), CCFs (lower right), bisec- 
tors (upper left), and bisectors velocity-span (lower left). The 
CCFs are clearly variable; the bisectors velocity-span diagram is 
composite: part of the data are spread horizontally over a large 
velocity range, and part are spread vertically, over a large range 
of span . The points that give the vertical bisectors velocity-span 
are those associated to the nightly high frequency RV variations; 
their bisectors are strongly variable in shape. The points with the 
low RV are associated to bisectors spans that are clearly shifted 
from the ones corresponding to higher velocities; this produces 
a shifted bisectors velocity-span. 
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Fig. 5. Simulation of a composite bisectors velocity-span dia- 
gram produced when adding a 100-Mj up companion on a cir- 
cular orbit, with a 120-day period around a 1.8-Mq pulsating 
star, HD 159492 (A7V; vsin/ = eOkms 1 ). The initial RV and 
bisectors velocity-span data are showed on the left; we see in 
particular high frequency (nightly) RV variations and bisector 
velocity-spans spread vertically; the simulated data are showed 
on the right. The bisectors velocity-span diagram on the right is 
clearly composite: both flat over a wide range of RV + vertical 
over a wide range of span values, similarly to HD 19545. 



called "anti-correlation"). For these objects, the ratio of the bi- 
sector span amplitude to the RV amplitude is found to be in the 
range 1-3. Figure [TJprovides an example of a star showing clear 
signatures of spots on the basis of the bisectors velocity-span 
diagram (HD 25457; F5V; vsin; = 25kms~ 1 ). As another ex- 
ample, the very nice case of HD 138763 can also be found in 
IDesort et al. 20071 

In the case of pulsations, the bisectors velocity-span values 
are spread over a much larger range than the RV and their vari- 
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Fig. 6. Example of a star whose RV variations are most probably 
due to binarity, HD 200761 (Al V; v sin i = 80 km s" 1 ). RV curve. 
This is a case where no CCF could be computed, and the binarity 
classification relies solely upon the RV curve. 
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Fig. 7. Example of a star with variable RV due to the presence 
of spots. HD 25457 (F5V; vsin / = 25kms~ 1 ). RV curve (upper 
right), CCFs (lower right), bisectors (upper left), and bisectors 
velocity-span (lower left). The bisectors velocity-span variations 
are clearly correlated to the RV variations. 
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Fig. 8. Example of a star with variable RV due to pulsa- 
tions. HD 159492 (A7V; vsin/ = 60kms _I ). RV curve (upper 
right), CCFs (lower right), bisectors (upper left), and bisectors 
velocity-span (lower left). The CCF and bisector variations are 
due to line deformations; the bisectors velocity-span variations 
are not correlated to the RV variations. 



ations are not correlated to the RV ones. In a (RV; bisectors 
velocity-span) diagram, the bisectors velocity-span values are 
spread vertically, and the ratio of the bisector span amplitude to 
the RV amplitude is large, typically > 3. Figure |8]provides an ex- 
ample of a pulsating star (HD 159492; A7V; v sin / = 60 km s~' ). 
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3. Results 

Given the variability criteria described above, 108 stars out of 
170 are found to be variable in RV, and 62 are found to be con- 
stant in RV within our precision limits. Table[T]provides relevent 
measurements on these targets: RV amplitudes and uncertainties, 
bisector velocity-span rms and uncertainties. 

3. 1 . Variability classification 
3.1.1. Stellar binaries 

20 stars are identified as binaries or candidate binaries with the 
criteria given in the previous section. More precisely: 

- 4 binaries are found on the basis of the x 2 criterium, namely 
HD 99453, HD 209819, HD2885 (Fig.S, HD 142629. 

- 6 stars show mostly flat bisectors velocity-span in a (RV; 
bisectors velocity-span) diagram: HD 11262, HD 68456, 
HD 41742, HD 116568, HD 216627, HD 54834. Their RV 
amplitude varies between 1 600 and 9 200 m s _1 . 

- 4 stars have composite, fiat+vertical bisectors velocity-spans 
in a (RV; bisectors velocity-span) diagram, together with a 
total RV amplitude dominated by the binarity effect. These 
pulsating binaries are: HD 220729, HD 12311, HD 19545 
(Fig.[4} and HD 1 12934. 

- Finally, 6 stars are classified as probable binaries on the 
sole basis of their RV variations: HD 158352, HD 177756, 
HD 158094, HD 2834, HD 200761 (Fig.|6]l, HD 1 16160. 

The stars are flagged in Table[T[ and an indication on the cri- 
teria which was used to identify them as binaries or possible bi- 
naries is also given. No attempt was made to further characterize 
the stellar companion once the binary status was established, and 
no more data were recorded on the objects. Their RV variations 
are given in Figure [9] 

Notes on some individual binaries (and potential binaries): 

- HD 11262 is associated to a ROSAT source by 
ISuchkov et al. 20031 



complicated than for pulsating stars. Our limited number of 
data does not permit to characterize the high-frequency pe- 
riod. 

HD 116160 was reported as an astrometric binary with ac- 



HD 54834: |Koen & Eyer 2002| (2002) reported this star as a 
photometric Hipparcos variable at a level of 0.0046 mag and 
with a frequency of 0.802 day -1 . Our data do not permit to 
confirm or infirm this frequency (not enough points, sam- 
pling not adapted). 

HD 68456 (Fig. [3]) was not reported as binary in the 
Hipparcos catalog from the photometric and astrometric 
points of view; it is classified by Adelman 2001 as one of 
the Hipparcos least variable stars. Goldin & Makarov 2007 
however provide an orbital solution to fit the Hipparcos 
astrometric data. The period found is 48 3 ±20 days, a ( ) = 



9 6 +26 



mas, eccentricity = 0. 12+"^, inclination = 131 + 16° 



cj = 103+^° and £2 = 171+^ 4 °. Fixing the period and eccen- 
tricity proposed by these authors, we tried to find a fit to our 
RV data. They happen to provide good fits assuming a mass 
of ^ lOOMj up for the companion. 

- HD 99453: |Baade & Kjeldsen 1997] questioned the previ- 
ously suggested SB2 status of this object on the basis of their 
data; we do confirm the SB2 status for this star. 

- HD 1 129 34 (A9V; v sin ; = 70 km s" 1 ): using Hipparcos pho- 
tometry, Handl er - 19991 reports this star as a new possible 
y Doradus candidate but with a "weak complicated signal", 
associated to a 0.8-day period. IdeCat et al. 20 06 did not find 
clear line-profile variations in their CORALIE data. From our 
data, the star is both pulsating and member of a binary sys- 
tem, which makes the line-profile variations indeed more 



celerating proper motion by Makarov & Kaplan 2005 



- HD 116568 was classified as one of the least variable stars 
with Hipparcos by lAdelman 20"0T1 |Baade & Kjeldsen 1997] 
report no variations in their ±0.5 km s" 1 spectroscopic sur- 
vey. The present data show that this star is a binary with an 
amplitude of at least 2750ms~'. It is also reported as an 
unresolved Hipparcos problem star by Masson 1999 and as- 
sociated to a ROSAT source by Suchk o~et al. 20031 

- HD 142629 is an astrometric Hipparcos binary. It was also 
recently reported for the first time as a spectroscopic binary 
bv lAntonello et al. 20"06l 

- HD 158352 was classified as a possible Herbig AeBe star by 
IThe et al. 19941 |Corporon & Lagrange" 1999 in a survey of 
RV variations among Herbig AeBe stars did not find varia- 
tions to a 5-10kms _1 level. This star was reported as being 
surrounded by a dusty disk by |Oudmaijer et al. 1992 and 
IMoor et al. 20061 give an age of 750±150Myrs for the sys- 
tem. 

- HD 177756 was classified as a possible ABootis 
star, as well as a possible SB (Farragiana et al. 2004 
IGerbaldi et al. 20031 . It is reported as one of the Hipparcos 
least variable stars (Adelman 200J}. 

- HD 200761 was reported as one of the Hipparcos least vari- 
able stars (Adelman 2001). 

- HD 209819 was also reported as one of the Hipparcos least 
variable stars (Adel man 200 II) . 

- HD 220729 is associated to a ROSAT source by 
ISuchkov et al. 20031 



3.1.2. Stars with planets 

One star, HD60532 (F6IV-V; B — V - 0.52) clearly reveals 
at the same time low-amplitude RV variations and flat bisec- 
tors velocity-spans diagram, indicative of the presence of two 
Jupiter mass companions with a high-confidence level. This 
star and the results of the fits of the RV curve is presented in 
IDesort et al. 200 8a Interestingly in the frame of the present pa- 
per, the periods of the detected planets are long (> 100 days). 
Hence, we get at least 1 % of F stars with long-period planets 
in our sample. This is much smaller than the predicted rate of 
^ 10 % for 1.5 M© stars by |Kennedy & Kenyon 2008 [ however, 
we are yet not sensitive to all range of masses and periods as will 
be showed in the last section. 



3.1 .3. Single stars: intrinsic variability 

We report in Table[T] the RV rms values obtained for each star, 
together with the associated uncertainties. Figure[l0]provides for 
all the stars except those identified as binaries the measured RV 
rms as a function of their B - V, the ratio RV rms/uncertainties 
(E/I) as a function of their B - V as well, and the (B - V; v sin i) 
diagram for the same objects. In the plots we have distinguished 
the 88 stars that are found to be variable according to the criteria 
defined above and those 62 found to be constant according to the 
same criteria. 

The E/I ratio varies between 1 .5 and a few tens; it is relatively 
smaller for stars with small B - V than for those with larger B - V . 
More quantitatively, the median value for this ratio computed 
on variable stars is 2.7 (resp. 4.4 and 5.2) for stars with B - 
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Fig. 9. Radial velocity curves of the identified or strong candidates or probable binaries (see text). 



V < 0.2 (resp. 0.2 < B - V< 0.4 and B - V> 0.4). Hence we 
detect more variable stars among stars with large B - V than 
stars with smaller B - V. We see moreover that the uncertainties 
generally increase with decreasing B — V. These results are not 
surprising and illustrate the fact that it is more difficult to identify 
variable stars when they have large uncertainties. In the frame of 
this study, it is important to keep in mind that our ability to detect 
variability generally decreases with decreasing B — V. 

Note also that the uncertainties increase with increasing 
vsinz'; we could actually verify that the uncertainties vary as 
vsku with a (vsku') a law where a = 1.5 ± 0.1, as predicted 
in lGalland et al. 2005al 

The percentage of variable stars depends on B-V in the fol- 
lowing way: 

- Most (85 %) of the 58 stars with B-V larger than 
0.4 are found to be variable and the RV uncertainty is 
2ms _1 (median value). Evenmore, 90 % of the 46 stars with 



B-V larger than 0.45, i.e., well beyond the instabil- 
ity strip, are found to be variable and their uncertainty is 
1.4 ms -1 (median value). We conclude then that at a level of 
precision of 2 m s~'or less, most of the stars with B-V larger 
than 0.4 are RV variable. 

- Among the stars with B-V between 0.2 and 0.4, the number 
of variable stars found is small, but this is due to a selection 
effect as known £Scuti and yDoradus stars were removed 
from our sample (see above). 

- Only 36 % of the 73 stars with B-V smaller than 0.2 are 
found to be variable. The percentage of variable decreases 
to 20% if we consider the 40 stars with B-V smaller 
than 0.1. For those stars with B-V between 0.1 and 0.2, 
we get as many variable as constant stars. The number of 
stars found to be constant according to our criteria increases 
then with decreasing B-V. However, we have seen that 
our ability to detect variable stars decreases with decreas- 
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Fig. 10. Top: RV rms measured for all stars but binaries with 
more than 6 spectra available as a function of B - V. Middle: 
Ratio RV rms/uncertainty for the same stars. Bottom: (B - V; 
v sin z) diagram for the same stars. Losanges indicate RV variable 
stars and square indicate RV constant stars. 



ing B — V. More quantitatively, the median uncertainty in 
the case of "constant" stars is 290 ms -1 whereas the me- 
dian uncertainty in the case of stars found to be variable is 
80ms -1 . Furthermore, the median uncertainty of constant 
stars is comparable to the median value of the standard de- 
viation of variable stars (265 ms -1 ). Hence we may stipulate 
that in fact most of the stars with B - V < 0.2 are probably 
RV variable. 

When the CCF and bisectors velocity-span criteria apply (in 
fact, whenever the bisectors velocity-span can be measured with 
a good or acceptable quality), we may try to further characterize 
the stellar variability. We find that in such cases, as expected, 
most of the variable stars with B — V smaller than 0.3 show signs 
of pulsations, whereas most of the variable stars with higher B - 
V show signs of spots. 

3.2. Variability of stellar origin and impact on planet 
detectability 

The "uncorrected" jitters, as given directly by the measured RV 
rms are provided in Table[T] for each star, together with the as- 
sociated uncertainties. Note that we prefer not to use the jitters 
corrected from the uncertainties, as sometimes done, as our main 
aim is to evaluate the impact on planet detectabiliity rather than 



to make stellar studies. Table|2] gives the computed median "un- 
corrected" jitters per bins of B - V. 

In lGalland et al. 2005al we had shown that the detection limit 
strongly depends on the star ST and its projected rotational- 
velocity; more precisely, the detection limit increases with ear- 
lier ST and/or larger v sin i. Thanks to the present data, we can 
in addition address the question of the impact of the stellar jitter. 

To estimate the detectable masses for a given star and a given 
period, we assume that a planet, supposedly on a circular orbit is 
detectable if the amplitude (2xK) of RV variations that it would 
produce is larger than 3x RV rms, where RV rms is the "uncor- 
rected" jitter actually measured. We will come back later on the 
validation on this assumption. We give for each star in Table 3 
the computed detectable limits assuming 3-day, as well as 10- 
day and 100-day periods. Figure 11 shows the detection limits 
for all stars for a 3-day period. For comparison, we also give 
in this figure the mass of the planet that would be detectable if 
the star is not active/pulsating (hence has no jitter) and the limit 
would then be set by the uncertainty (hypothesis 2K = 3x un- 
certainty). The ratio of the two values is E/I. As previously seen, 
this ratio is larger than 1 .5 and may be quite high; the impact of 
the jitter on the detectable masses is therefore non negligeable. 

Several comments can be made: 

- the achievable limits fall into the planetary domain for a large 
number of stars: more precisely, in 137 out of 150 stars, i.e., 
91%, the detection limit for a 3 -day period falls within the 
planetary domain. For the remaining stars the limit falls well 
into the BD domain with masses up to 54Mj up . When con- 
sidering a 10-day (resp. 100-day) period, we find that we can 
reach the planetary mass domain for 124 stars, hence 83 % 
(resp. 92 stars, hence 61 %). For a 10-day period, the limit 
for all remaining stars but one fall into the BD regime; for 
the 100-day period, the limit for all remaining stars but 6 fall 
in the BD domain. 

- as expected, the median of the detection limits generally im- 
prove with increasing B - V, from 10Mj up for B - V be- 
tween -0.1 and 0, to 5 M Jup for B - V between and 0.3, to 
0.05 M Jup for B - V between 0.5 and 0.6 (for a 3 day period) 
(see Table|2]i. Noticeably, for stars with B - V > 0.3, individ- 
ual detection limits may be as low as 0.02 Mj up and for stars 
with B — V <, 0.3, individual detection limits may be as low 
as 0.5 Mj up . For a ten day period, these numbers become re- 
spectively: 15, 7 and 0.08 M Jup ; for a 100 day period, 31, 16 
and 0.17M Jup . Also, noticeably, the detection limits steeply 
improves at B - V = 0.3. 

- the "uncorrected" jitter varies a lot from one object to the 
other. Therefore the general conclusion that the detection 
limits improves with increasing B - V may not apply when 
considering individual objects: for instance, the two stars 
HD 50445 (A3V; B-V = 0.18) andHD 63847 (A9V; B-V = 
0.3) have similar projected rotational-velocities (vsinz - 
90kms _1 ) and very different levels of activity, with an RV 
rms of 66 ms -1 and 794 ms -1 respectively. When we take 
this "uncorrected" jitter into account, the detection limit is 
lMj U p (P = 3 days) and 1.5M Jup (P = 10 days) around the 
A3V star whereas the detection limit is about 10 times higher 
for the A9V star. 

We conclude then that planets can indeed be found around a 
wide range of stars with B-V larger than -0. 1, even taking into 
account their jitter. The achievable detection limit of such early 
type stars cannot be predicted given only the star properties (ST, 
v sin z), but requires to record data to estimate their level of jitter. 
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Table 2. Median values of RV rms (second row) and RV uncertainties (third row) for stars with more than 6 spectra (3 epochs) 
available. The number of stars per bin of B - V is also given (first row). Rows 7-9 give the median per bin of B - V of the achievable 
detection limits deduced from the measured rms for each star, expressed in Jupiter mass. Three periods are considered: 3, 10 and 
100 days (see Sect.|3]l. Note that we assume that the planet, supposedly on a circular orbit, is detectable if the amplitude of the RV 
variations is larger than 3 x RVrms. Rows 7-9 give the percentage of stars for which the detection limit, given the measured RV 
rms, fall in the planetary domain. 
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Fig. 11. Achievable detection limits for all stars but binaries tak- 
ing into account their actual "uncorrected" jitter (losanges) or 
the measured uncertainties (squares), and assuming a planet on 
a circular 3 -day orbit. 

Note: we note that of course, the measured "uncorrected" 
jitter provides a reliable limit to planet detection only when this 
jitter is due to stellar activity in general, and not to companions. 
Would a companion be present, its contribution to the RV vari- 
ability would have to be removed in order to estimate the impact 
of the stellar activity. 

4. Planet detection limits of the present survey 

4.1. Estimation of the achieved detection limits 

We try now to estimate the detection limits reached by the 
present survey, taking into account the actual RV curve. For each 
star, we then compute the detection limit (companion mass) as a 
function of its period. To do so, we consider a planet with a given 
mass and with a given period (the orbit is assumed to be circu- 
lar). For any couple (Mass; Period) we generate a large number 
of keplerian orbits, assuming different times of passage at peri- 
astron (To). For each orbit, we compute the expected radial ve- 
locities at the times of the actual observations. We add a noise 
(random value between t RV uncertainty), where RV uncertainty 
is the uncertainty measured on the RV data. We then get a virtual 



set of RVs, which takes into account the star properties (in par- 
ticular, its ST and rotational velocity, through the uncertainties 
and SN). We compute then the standard deviation of the virtual 
RVs points. For a given (Mass; Period), the distribution of all the 
standard deviations (corresponding to differents 7q) obtained is 
gaussian. We compute then the average value of the distribution 
of the virtual standard deviations. We consider that a planet with 
a given (Mass; Period) is detectable if the standard deviation of 
the real RV values is smaller than the average value of the virtual 
standard deviations. We determine the confidence level (or de- 
tection probability) associated to such an orbit by comparing the 
standard deviation of the virtual distribution with the difference 
between the standard deviation of the real RV measurements and 
the average value of the virtual standard deviations. 

In practice, for a given object, we explore 200 periods in 
the range 0.5 to 1000 days, and 100 planet masses in the range 
(Mmin; 100Mj up ) where Mmin corresponds to the achievable 
mass given the measured uncertainty. For a given (Mass; Period), 
we explore 1000 Tq. We checked that increasing the number of 
periods, or planet masses and or initial To does not significantly 
impact the results. 

For each (Mass; Period) couple, we obtain thus a detection 
probability. In a (Mass; Period) diagram, we can then identify 
the domain where a planet with a given mass and period should 
be detectable if present, with a given level of confidence. This 
defines then a domain in which we can exclude the presence of 
a planet with a given level of confidence. We will consider two 
levels of confidence: lcr, (i.e., a 68.2 % probability) and 3cr (i.e., 
99.7 % probability). 

4.2. Sensitivity of the survey and first constrains on early 
type stars 

The sensitivity of our survey is a consequence of the number 
of data available and on the temporal sampling of the data. We 
kept only those stars (107 objects), found to be either constant or 
variable, for which we got more than 12 data points (6 epochs). 
Also, given the data at hand, we limited the range of periods 
investigated between 1 days and a few hundred days. We report 
in Table 3 the achieved limits (68.2 % and 99.7 % probabilities) 
obtained for each of the 107 stars considering three periods: 3, 
10 and 100 days. 

We give also in Fig.[l"2"]examples of achieved detection lim- 
its (68;2 % and 99.7 % probabilities) as estimated with the previ- 
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ously described simulations. We have also plotted on the figures 
the achievable detection limits taking into acount the jitter, as 
defined in the previous section, as well as the detection limits 
corresponding to the measured uncertainties. We recall that the 
last two cases (achievable limits) do not take into account the 
actual temporal sampling of the data, conversely to the detection 
limits computed with our virtual realisations. 

When enough data are available, the actual detection limits 
fall close to the achievable limits obtained assuming the 3 x RV 
rms threshold for the amplitude of RV variations, as can be seen 
in Fig. [12] This justifies the choice of the threshold adopted in 
the previous section to estimate the achievable limits. 

Note that we sometimes endup with high detection lim- 
its that fall outside the investigated range of masses, i.e., > 
100 Mj up when we consider a period of 100 days whereas the de- 
tection limits for a 3 or 10 days period are close to the achievable 
limits. This corresponds to cases where the temporal sampling is 
not adapted to explore such a long period (see for instance the 
case of HD 33256, in Fig. 12 1. It sometimes, but much less fre- 
quently also happens that the detection limit for a 10 days fall 
outside the investigated range of masses whereas the detection 
limits for a 3 days period is close to the achievable limit. In fact, 
our survey allows to search mostly for short-period planets (typ. 
a few days); 10-100 day periods are not always sampled enough 
to get interesting results (especially on early-type stars) and > 
100-day periods are not properly sampled to get interesting re- 
sults. We will therefore in the following discuss only periods < 
100 days. Finally, one has to note that in some cases, we get a 
99.7 % probability detection limit out of the investigated range, 
where as the 68.2 % limit fall well into the investigated range. 
This happens generally when the number of data is the lowest: 
12 or 14. 

To study the impact of B - V on the present results, we com- 
puted the percentage of stars for which the achieved detection 
limits (68.2 % and 99.7 % probabilities) fall in the planetary or 
BD domains per bin of B- V, considering a 3-day, a 10-day and a 
100-day periods. We also computed the median of the achieved 
detection limits (considering 68.2 % and 99.7 % probabilities) 
per bin of B- V for such periods. The results are given in Table|4] 
In order to allow comparison between the limits obtained with 
these two probabilities and with the achievable ones, we consid- 
ered for the computation of the median values only those stars 
for which both the 68.2 % and 99.7 % probabilities detection 
limits fall into the investigated range of masses. Finally, one has 
to note that for the earliest type stars, the number of objects per 
bin is quite small, so one has to be very cautious with the asso- 
ciated statistics. We can see that: 

- if we consider a 3-day period, the achieved limit at 1 (resp. 
3) cr fall in the planetary domain for 90 % (resp. 81 %) of the 
stars. This percentage is quite comparable to the one found in 
Sect. 3. It increases from 75 (resp. 25) % for the earliest type 
stars to 100 (resp. 100) % for stars with B - V larger than 
0.3. Also, the median of the achieved limits at 1 (resp. 3) cr 
decreases from 7 (resp. 7) M Jup for the earliest-type stars to 
0.08 (resp. 0.3) M Jup for the latest-type stars. Moreover, the 
steep step seen in Sect.|3]in the detectable masses at B - V — 
0.3 is also clear. 

- if we consider a 10-day period, the achieved limit at 1 <x 
(resp. 3) fall in the planetary domain for 82 % (resp. 67 %) of 
the stars. This percentage increases from 50 (resp. 25) % for 
the earliest type stars to 100 (resp. 100) % for stars with B - 
V larger than 0.4, with however an exception on the [0.2; 0.3] 
range where it decreases back to 50 %. Also, the median of 



Table 5. Detection limit for 50 % and 90 % for different peri- 
ods. The last column give the number of stars considered to esti- 
mate these detection limits. Note that only stars with more than 
12 measurements (6 epochs) were considered. Binaries were ex- 
cluded. Also note that numbers outside the planetary mass do- 
main (> 13 Mj up ) are not given. 
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the achieved limits at 1 (resp. 3) cr decreases from 12 (resp. 
24) Mj up for the earliest-type stars to 0.1 (resp. 0.7) for the 
latest-type stars, with however an exception in the range [0.2; 
0.3] range as regards the 99.7 % probability. Again the steep 
step is observed at B - V = 0.3. 
- if we finally consider a 100-day period, the achieved limit at 
1 (resp. 3) cr fall in the planetary domain for 54 (resp. 35) % 
of the stars. We note that this percentage is smaller than the 
one obtained in Sect. 3 and attribute the discrepancy to the 
actual temporal sampling and the relatively small amount of 
targets investigated yet. The median of the achieved limits at 
1 (resp. 3) cr decreases from 19 (resp. 34) Mj up for stars with 
B - V > 0.0 to 1. (resp. 1.9) M Jup for the latest-type stars. 
Again the steep step is observed at B - V — 0.3. 

Finally, we give the probability not to detect planets of a 
given mass and with a given period (3, 10, 200, 500 days) around 
stars with a given spectral type and v sin i. The results are sum- 
marized in Table[5] We see that as expected, for a given proba- 
bility, the limits globally decrease with increasing B-V and de- 
creasing v sin i. 

Obviously, the statistics provided by our survey is still poor 
on early type stars, and still very limited on the latest type stars. 
Concerning the latter, we note that if we consider the 41 objects 
with B - V > 0.4, i.e., well beyond the instability strip, we find 
that less than 24 % of stars host planets with masses equal to 
0.5 Mj up or more; less than 5 % host planets with masses equal 
to 1 Mj up or more on a 3-day period. For a 10-day period, we 
find that less than 42 % host planets with masses > 0.5 Mj up , and 
less than 20 %) host planets with masses > 1 Mj up . The compari- 
son between achieved and achievable detection limits shows that 



A.-M. Lagrange et al.: Extrasolar planets and brown dwarfs around A-F type stars 



11 



100.0 i—^ 



10.0 



HD30652 - F6V - 16 km s" 



0.01 




i 1 100.0 



10.0 



1.0 



10.0 



100.0 



HD32977 - A5V - 100 km s" 



0.01 
1000. 1.0 




10.0 



100.0 



1000. 



HD33256 - F2V - 10.0 km s" 
100.0 ir 



10.0 - 



3 



„ o.oi 




HD38393 - F7V - 8.4 km s" 
100.0 r— <~i 1 r 



10.0 - 




0.01 



1.0 



10.0 



100.0 



1000. 



1.0 



10.0 



100.0 



1000. 



B 100.0 



o 
S 

Q 



HD50445 - A3V - 91 km s" 




0.01 



1.0 



10.0 



100.0 



100.0 



10.0 



1.0 



0.1 



HD63847 - A9V - 94 km s" 




x * x 



-x 



0.01 
1000. 1.0 



I ,_i i L 



J I i L 



10.0 



100.0 



1000. 



100.0 



10.0 



HD1 15892 - A2V - 90 km s" 




0.01 



100.0 



HD147449 - FOV - 83 km s" 



1.0 



10.0 



100.0 



0.01 
1000. 1.0 




10.0 



100.0 



1000. 



Period [days] 

Fig. 12. Detection limits. X-axis: periods (days). Y-axis: Detection limit (M/Mj up ). Curve: detection limits actually achieved in the 
present survey; plain curves correspond to 99.7 % detection probability, and broken curves to 68.2 % probability. Note that the RV 
data were beforehand averaged over one day. Straight line: achievable detection limits assuming a 3x rms threshold for the planet 
amplitude. Line with crosses: achievable detection limits assuming a 3x uncertainty threshold for the planet amplitude. 
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Table 4. Percentage of stars for which the achieved detection P=68.2% or P=99.7% limits fall in the planetary or BD/planet domains. 
Three periods are considered: 3, 10 and 100 days. Median values of the achieved detection limits, expressed in Jupiter mass, per 
bin of B - V. For the these median values, we took into acount only the stars for which both the P=68.2% and P=99.7% limits fall 
into the investigated companion mass range. The corresponding number of stars is given as well as the median achievable limits for 
those stars. 
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P = 3 days 
Number of stars 
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17 
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22 
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Achieved detection limit (P=68.2%) 


6.9 


5.3 


5.1 


3.2 


0.5 


0.25 
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Achieved detection limit (P=99.7%) 


6.9 


12.1 


7.3 


6.3 


1.0 


0.3 


0.3 


Achievable detection limit 


6.9 


5.5 


4.2 


2.5 


0.4 


0.2 


0.06 


P = 10 days 
Number of stars 


4 


14 


20 


7 


11 


22 


17 


Achieved detection limit (P=68.2%) 


11.7 


8.0 


7.2 


4.4 


0.9 


0.4 


0.1 


Achieved detection limit (P=99.7%) 


24.0 


11.8 


10.2 


17.4 


2.0 


0.8 


0.7 


Achievable detection limit 


10.3 


6.9 


5.6 


4.1 


0.6 


0.3 


0.1 


P = 100 days 
Number of stars 


1 
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8 


13 
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Achieved detection limit (P=68.2%) 
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17.3 


14.2 


5.2 


1.1 


1.0 


Achieved detection limit (P=99.7%) 


(17.6) 


34.4 


32.9 


33.8 


20.4 


2.4 
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(8.5) 


15.0 


10.5 


8.9 
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there is still room to significantly improve those statistics (thanks 
to new data points). 

The present statistics certainly does not allow quantitative 
comparisons with late type dwarfs, which have been surveyed 
by several groups since more than 10 years, or with giant or sub 
giant stars, because in that case of the lack of data for both mas- 
sive dwarfs and (sub-)giants. 

Concerning the presence or absence of hot Jupiters around 
massive stars, we note that the planets found so far in our sur- 
vey are located at about 0.7 AU or more from a 1.4 Mq star; 
this separation corresponds to that of the closest planet found 
around giant stars. We also detected recently a planet orbiting at 
0.6 AU from a dwarf with a similar mass in the northern hemi- 
sphere dDesort et al. 2008bb . Due to the still limited amount of 
data available, it should not be concluded however that there 
are no planets closer to massive dwarfs. We also recall that a 
few short period planets have been found around 1 .4 M0 stars 
through transits. 

5. Conclusion 

Based on the observation of a large number of A-F type stars 
(170), we have been able to measure their jitters, and derive for 
the first time estimations of the detection limits that can be ex- 
pected on average on those stars with B - V in that range [-0.1; 



0.6] (once previously known £ Scuti and yDoradus stars are re- 
moved) for 3 periods: 3, 10 and 100-day periods. We have shown 
that at the precision provided by HARPS, most of the stars are 
variable in RV, and the impact of the RV jitter, due to either spots 
or pulsations is generally non negligeable on planet detectability. 
However, assuming that planets are detectable if the amplitude 
of the induced RV variation is larger than 3xrms (this threshold 
defines the achievable detection limits, which depends on the star 
and the spectrograph used), we have shown that even when tak- 
ing into account the jitter, giant planets can still be found around 
these stars in most cases. This is not only true for the stars with 
B - V > 0.3 for which we can find either short or long period 
planets, with masses as small as 0.02 M Jup (case of short period) 
for the latest type stars but also for dwarfs with B - V < 0.3: for 
such stars short-period planets can still be found around those 
with relatively low projected rotational-velocity and low level of 
activity, with masses down to 0.5 Mj up (best case). This survey 
has allowed to identify for the first time the stars that are best 
suited for further searches for planets around massive dwarfs. 

We have also shown that given the data available, the present 
survey is sensitive to short-period planets (hot Jupiters) and only 
partially sensitive to larger periods (up to 100 days). We found in 
particular one 2-planet system with periods larger than 100 days 
around one late-type star. Whenever possible (107 stars), we 
computed for each star the detection limits actually achieved and 
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showed that when enough data are available, the achieved detec- 
tion limit is set by the 3 xrms threshold. We indeed reached such 
limits for early-type as well as for late-type stars. We finally de- 
rive first estimates of the presence of short-period planets around 
these A-F stars. We showed for instance that less than 5 % of 
the latest-type stars (B — V > 0.4) host P=3 days-period planets 
with masses 1 Mj up or more. Such statistics are not constraining 
enough to allow interesting comparisons with later type stars or 
with model predictions; but as soon as more data become avail- 
able, the statistics can be straightforwardy improved. 

Finally, we note that to compute these detection limits, we 
did not try to average out the spectra over timescales associated 
to the frequencies of intrinsic stellar variations. This approach 
would of course allow to decrease significantly the detectable 
masses. As it would require lots of telescope time, it should be 
probably kept for stars with the highest scientific interest. 
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Table 1. Stars properties and measurements. "RV rms" (resp. "span rms") stands for the rms of the measured radial velocities (resp. bisector 
velocity-spans); "RV amp" (resp. "Span amp.") stands for amplitude of the measured radial velocities (resp. bisector velocity-spans); "RV una" 
(resp. "Span una") stands for the average uncertainties associated to the RV (resp. bisector velocity-spans) data. Bisector flags: G: good quality; 
B: bad quality; A: acceptable quality. Binary types: X refers to stars identified as binaries based on a high^ 2 (> 10); B refers to binaries identified 
via a flat or a composite bisector, and V refers to stars regarded as binary candidates, based on the sole amplitude of their RV variations (see text). 
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Table 3. Detection limits, either achievable or achieved with a 68.2% or a 99.7% probability in the present survey (see text). Note that only those 
detection limits less than 0.05 Jupiter mass are given with 2 digits. 
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